Buckling of Carbon Nanotube-Reinforced Polymer Laminated Composite Materials 
Subjected to Axial Compression and Shear Loadings 

J. C. Riddick 1 

Vehicle Technology Directorate, U.S. Army Research Laboratory, Hampton, VA 23681 

T. S. Gates 2 

NASA Langley Research Center, Hampton, VA 23681 
S.-J. V. Frankland 3 

National Institute of Aerospace, 100 Exploration Way, Hampton, VA 23666 


ABSTRACT 

A multi-scale method to predict the stiffness and stability properties of carbon 
nanotube-reinforced laminates has been developed. This method is used in the prediction 
of the buckling behavior of laminated carbon nanotube-polyethylene composites formed 
by stacking layers of carbon nanotube-reinforced polymer with the nanotube alignment 
axes of each layer oriented in different directions. Linking of intrinsic, nanoscale-material 
definitions to finite scale-structural properties is achieved via a hierarchical approach in 
which the elastic properties of the reinforced layers are predicted by an equivalent- 
continuum modeling technique. Solutions for infinitely long symmetrically laminated 
nano tube-reinforced laminates with simply-supported or clamped edges subjected to axial 
compression and shear loadings are presented. The study focuses on the influence of 
nanotube volume fraction, length, orientation, and functionalization on finite-scale 
laminate response. Results indicate that for the selected laminate configurations 
considered in this study, angle-ply laminates composed of aligned, non-functionalized 
carbon nanotube-reinforced lamina exhibit the greatest buckling resistance with 1% 
nanotube volume fraction of 450 mn uniformly-distributed carbon nanotubes. In addition, 
hybrid laminates were considered by varying either the volume fraction or nanotube 
length through-the-thickness of a quasi-isotropic laminate. The ratio of buckling load-to- 
nanotube weight percent for the hybrid laminates considered indicate the potential for 
increasing the buckling efficiency of nanotube-reinforced laminates by optimizing 
nanotube size and propotion with respect to laminate configuration. 


INTRODUCTION 

Because of many unique mechanical and electrical properties, carbon nanotubes 
(CNT) have been recognized as particularly well suited for the development of 
nanocomposite structural materials. For applications that require thin gauge components, 
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load carrying capability, and utilize minimum structural weight as a key design objective, 
one choice is to incorporate the CNT in a thin, polymer film or layer, to form a 
nanocomposite material. Potential applications for such nanostructured films include 
surface coatings, MEMS sensors, actuators, and gossamer or ultra-light aircraft. 
Experimental characterization using thermogravimetric analysis, mechanical tests, and 
scanning electron or probe microscopy indicates that the addition of CNT enhances 
mechanical properties of polymer films with respect to hardness, stiffness, and flexural 
strength. 

Recent advances have lead to the synthesis of ultra-thin layered nancomposites. A 
major purpose of layering, or lamination, is to tailor the directional dependence of 
strength and stiffness of a material to match the loading environment of a structural 
element. Reference 1 discusses experimental issues related to fabrication of 
nanolaminates such as adsorption and dispersion. In contrast to ultra-thin carbon-fiber 
composites, CNT-polymer composites are considered nanostructured in the sense that the 
primary constituents have structure that can be readily described at the nano-scale. These 
material descriptions then provide the basis for the development of the underlying 
structure-property relationships that can be subsequently used in analysis models to tailor 
or enhance performance of structural components. To be truly useful, the structure- 
property relationships and related models must span multiple length scales and be 
versatile enough to be used in parametric studies that influence material development and 
design. 

Buckling resistance is an important consideration in the design of thin-walled, 
lightweight aerospace vehicles and structural components. Reducing structural weight 
while at the same time increasing buckling capacity make composite materials, which 
exhibit high strength-to-weight and stiffness-to-weight ratios, ideal for aerospace vehicle 
structural components. Nanotube-reinforced laminated composites represent a new class 
of structural materials to be exploited for potential development of a variety of new, 
advanced applications. Ideally, longer nanotubes in high volume concentrations will offer 
the greatest mechanical advantage in enhancing bulk constitutive properties of matrix 
material. However, the current state-of-the-art in manufacturing CNT-polymer 
composites suggest that cost and processing issues will drive the future development 
CNT-reinforced laminated composite materials. 

Structural tailoring of laminated composite plates has been studied extensively. 
Nemeth in reference 2 provides an exhaustive parametric treatment of the buckling 
response of symmetrically laminated fiber-reinforced panels. The detailed studies 
contained in this reference culminates in a set of design curves based upon the buckling 
behavior of infinitely long laminates. Buckling results for infintely long plates are 
important to designers because they provide the lower bounds to the festoon buckling 
curves of finite-length rectangular plates that are particularly useful in preliminary design 
studies. The nondimensional buckling parameters and equations presented by Nemeth 
have been used by Odegard, et. al., 3 as part of a multi-scale analysis based upon a 
hierarchical approach used to determine the buckling resistance of single layer 
ortho tropic composite plates subjected to uniform axial compression, uniform shear, or 
pure in-plane bending for CNT-reinforced composite plates. In particular, plates made of 
a polyethylene matrix material reinforced with single-walled carbon nanotubes were 
analyzed. The multi-scale analysis was used to quantify the relative influence of 



molecular-level structural items defined as nanotube length, volume fraction, orientation, 
and functionalization on the buckling resistance of infinitely long plates. Each of the four 
molecular attributes considered may be controlled at material synthesis. The study found 
that randomly-oriented, non-functionalized nanotube-reinforcement would be 
recommended for single-layer plates when the design criteria is focused on providing the 
largest possible value for in-plane buckling resistance. 

In this paper, a multi-scale analysis will be employed to predict the stiffness and 
stability properties of CNT-reinforced laminated composite materials. The objective of 
the present study is to extend the existing multi-scale analysis method of Odegard, et. al., 3 
to predict the stiffness and stability properties of a CNT-reinforced laminated composites. 
The multi-scale analysis concept li nk s nanomechanics results from molecular dynamics 
(MD) simulations to the Mori-Tanaka micromechanics model via an equivalent- 
continuum modeling approach in order to predict bulk elastic properties for polymer 
lamina reinforced using single-walled carbon nanotubes (CNT). Herein, the CNT- 
dependent properties derived from the equivalent-continuum modeling approach are 
incorporated into traditional Classical Lamination Theory (CLT) to obtain a first-order 
prediction of macromechanical stiffness properties for a given CNT-reinforced laminate 
configuration. This method is used to facilitate the prediction of the buckling behavior of 
infinitely long symmetrically laminated CNT-reinforced polyethylene composites to be 
used in preliminary design studies. The symmetrically laminated CNT-reinforced 
polymer composites are formed by stacking layers with the principal material axes of 
each layer oriented in different directions, such that each layer located above the laminate 
midplane has a corresponding layer below the midplane, with the same lamina properties. 
For the present study lamina material properties are defined as thickness, nano tube 
volume fraction, length, alignment, and orientation. 

Nondimensional parameters and equations that define the stiffness and stability 
properties of laminated plates are applied to selected laminates subjected to uniform axial 
compression and shear loadings. Buckling solutions for infinitely long laminates with 
simply-supported or clamped edges are presented. Finally, hybrid laminates are formed 
by varying molecular properties through-the-thickness of the laminate. The buckling 
response of selected hybrid laminates is presented in terms of a buckling-to-weight 
percent nano tube ratio. Here, buckling load is normalized by nano tube weight percent in 
an effort to assess the potential for tailoring the laminate response in a cost effective 
manner with respect to laminate configuration. 


PROBLEM DESCRIPTION 

A description of the multi-scale analysis is presented here. First, the material system 
to be analyzed is described. Then, analysis details at each length scale are outlined. 
Finally, the details of the hybrid laminate definition and nomenclature are presented. 

Materials 


The material system used in this study was a carbon nanotube-reinforced polymer 
composite. The carbon nanotube material was modeled as a single- walled (10,10) 



nano tube of radius 6.78 A and variable length. The polymer matrix material was assumed 
to be isotropic, amorphous polyethylene matrix with a representative Young’s modulus 
and Poisson’s ratio of 0.9 GPa and 0.3, respectively. 

Nanomechanics 

The nano tube-reinforced CNT-polymer laminates considered in the present study are 
considered nanostructured because the primary constituents have structure that must be 
analyzed at nano-scale. A nano-scale analysis 3 is employed to develop underlying 
structure-property relationships that can be used in analysis models to tailor or enhance 
performance of structural components. Highlights of the nanomechanical analysis are 
presented in this section. 

Previous work 3 defined four nanoscale structural items as the primary molecular 
attributes, controlled at material synthesis, which could influence final properties of the 
composite. These four attributes are nanotube length, L, nanotube volume fraction, v/, 
nanotube orientation, 0, and nanotube functionalization. Nanotube functionalization is 
defined as the formation of a chemical covalent bond between the nanotube and polymer 
directly. The hierarchical approach employed in the present study allows the influence of 
these four intrinsic attributes on bulk constitutive properties to be quantified. 

The hierarchical modeling method relies on three major steps. First, molecular 
dynamics (MD) simulations are conducted using a suitable representative volume 
element (RVE) of the nanoscale molecular structure of the CNT-polymer system. 
Second, an equivalent-truss model is developed, using finite-element analysis, to link the 
MD results to an equivalent-continuum model. Finally, the equivalent-continuum model 
is developed to model the nanotube as an effective fiber. The effective mechanical 
properties of the molecular system are determined by equating the strain energies of the 
molecular and equivalent-continuum models. The equivalent-continuum-effective-fiber 
results yield the nine independent elastic constants used to describe the overall of the 
behavior of the effective fiber: three elastic axial stiffnesses, C//, C?/, and C?/; the 
plane-strain bulk moduli, K 2 /, Kii, and Kjf; and the three elastic shear stiffnesses, Cj, 
Cj/, and CrJ. Where, the subscripts 1, 2, 3 indicate principal material directions; 
subscripts 4, 5 and 6 are indices of the stiffness tensor written in Voigt’s notation; and the 
superscript/ denotes that these properties are associated with the effective fiber. 


Table 1. Equivalent-Continuum Properties (GPa). 



Non-functionalized 
effective fiber 

Functionalized 
effective fiber 

Cj/ 

548.7 

487.7 

C 2 1 

16.8 

24.5 

Csl 

16.5 

20.6 

k 2 / 

14.8 

19.5 

K,/ 

149.3 

137.1 

Kj 

149.2 

138.7 

cj 

7.1 

12.7 

cj 

144.0 

155.4 

cj 

144.9 

137.0 



The nine independent elastic properties given in Table 1 are determined by applying 
nine identical sets of boundary conditions to the equivalent-truss model and the 
equivalent-continuum model. The strain energies of the models are equated by adjusting 
the elastic properties in a parametric fashion. Table 1 presents the nine independent 
parameters, where functionalized and non-functionalized nanotube systems were 
modeled as effective fibers. 

Micromechanics 

Constitutive models governing the behavior of nanotube-reinforced polymer lamina 
were obtained from a micromechanical analysis by using the mechanical properties of 
effective fibers. 3 In this analysis, the polymer molecules near the polymer/nanotube 
interface were included in the effective liber, and it was assumed that the polymer matrix 
surrounding the effective fiber had mechanical properties equal to those of bulk 
polyethylene. Perfect bonding between the nanotube/polymer effective libers and the 
surrounding polymer matrix was also assumed. 

A micromechanics-based Mori-Tanaka method was used to predict the elastic 
mechanical properties of the composite material. 3 For this method, the overall elastic- 
stiffness tensor of the nanotube-reinforced lamina composed of CNT embedded in a 
matrix material is 

c=(' , m c m +v / ^cyr^))(v m i+Vf( , iy)) -1 (1) 

where v/ and v m are the effective fiber and matrix volume fractions, respectively, C) and 
C m are the stiffness tensors of the effective fiber and the matrix respectively, I is the 
identity tensor, the angle brackets indicate an effective-fiber orientation average, and Tf is 
the dilute strain-concentration tensor of the effective fiber, and is given by 



where S is the Eshelbyl3 tensor. 

Lamina properties, namely the Young’s and shear moduli are determined using the 
homogenized-composite stiffness tensor, C. The shear moduli of the composite material 
( G23 , G31, G12 ) are simply equal to the shear-stiffnesses (C44, C55, Cm), the Young’s 
moduli (Ei, E 2 , E3) were calculated by using the components of the compliance tensor of 
the nanotube-reinforced lamina detennined by inverting the composite stiffness tensor 
C . 15 For composite lamina composed of aligned CNT, the subscripts of the Young’s and 
shear moduli indicate the principal direction associated with the quantity. Where the 1- 
direction is parallel to the longitudinal axis of the nanotube, the orthogonal 2- and 3- 
directions are transverse to the logitudnal axis of the nanotube. For the random-oriented 
composite lamina, isotropic mechanical properties are completely described by the 
Young’s modulus, E (E = Ej = E 2 = E3 ) and shear modulus, G (G = G23 = G31 = G 72 ). 

For the present study, lamina properties were determined using the elastic stiffness 
components, volume fraction, length, and orientation of the effective fiber properties in 
Eq. (1). The calculations were performed for perfectly aligned and three-dimensional 



randomly oriented effective fibers. Symmetrically laminated CNT-reinforced polymer 
composites are formed by stacking lamina with the principal material axes of each layer 
oriented in different directions in a manner that maintains the symmetry about the 
laminate midplane, as depicted in Figure 1. 

Macromechanics 

The buckling analysis is based on the classical Rayleigh-Ritz variational method and 
is derived explicitly in terms of the nondimensional parameters defined in Reference 2. 
Deriving the analysis in this manner produces equations well-suited for parametric 
studies. The nondimensional parameters used are given by 



( 3 ) 

( 4 ) 

( 5 ) 


Where, the parameter /3 characterizes plate bending orthotropy, and 8 and y characterize 
plate bending anisotropy. The subscripted D terms appearing in the equations are the 
plate bending bending stiffnesses of Classical Lamination Theory (CLT) 4 defined as 

i) 

J k=i 

where k, N, z*, and z*./ are defined in Figure 1, Q, } represent the transformed reduced 

stiffnesses for each lamina as defined in Reference 4. For aligned nano tubes, the 
transformed reduced stiffness matrix is dependent on the nanotube orientation angle and 
the four independent elastic constants Ej, Ej, Gn, and Vn. For lamina reinforced by 
randomly-oriented nanotubes, the transformed reduced stiffnesses reduce to the isotropic 
formulation presented by Odegard, et. al. 3 Note that all of the parameters defined by 
equations (3) through (5) depend on plate bending stiffnesses which can be shown to be 
functions of the lamina elastic properties and nanotube orientation. The lamina 
constitutive properties and nanotube orientation are functions of the intrinsic molecular 
properties. Therefore, the bending stiffness, and furthermore the parameters of equations 
(3) through (5), may be expressed as functions of the intrinsic molecular properties. 

In addition to /?, 8, and y, two additional nondimensional quantities are used to 
characterize buckling resistance, and they are given by 



( 7 ) 


K x 


K V b 2 
''J D\\Di2 


K s = , ' ( 8 ) 
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The quantities K x and K s are referred to as the axial compression buckling coefficient and 
the shear buckling coefficient, respectively. Each of the loading conditions associated 
with these buckling coefficients is shown in Figure 2. The compressive and shear stress 
resultants corresponding to the loadings are denoted by N x and N xy , respectively. All 
stress resultants shown are positive-valued. The superscript cr denotes the critical value 
of the applied load associated with the onset of buckling. Finally, the laminate width, b, is 
defined in Figure 2. 

A set of simplified buckling equations are presented in Table 2 which were 
determined by performing a least squares fit to response curves determined by a special- 
purpose Rayleigh-Ritz analysis. Nemeth 2 provides the following simplified linear 
expressions for the buckling coefficients: 


Table 2. Buckling coefficients for in-plane loading. 



Simply-supported 

clamped edges 


edges 


K x 

2 + 2/3 

4.602+2.359 /3 

K 

3.629+ 1.644 /3 

6.493 + 2.414 /3 


The equations for the buckling coefficients given in Table 2 depend only on 
nondimensional parameter (3. These equations arise from an analysis that neglects the 
effect of laminate bending anisotropy. According to Reference 5, the bending anisotropy 
is expected to be negligible when the maximum value of nondimensional anisotropy 
parameters 6 and y are less than 0.2. A potential benefit of certain symmetric laminate 
configurations is a high degree of bending anisotropy. Such categories of laminates could 
potentially give rise to 6 and y values in excess of 0.2. Therefore, care must be taken in 
choosing laminates to be analyzed using these buckling equations. 

The potential tailoribility of the buckling behavior of nanotube-reinforced 
polymer composites will be assessed by studying the effects of fiber orientation, stacking 
sequence and intrinsic molecular properties on the nondimensional parameters and 
buckling coefficients. In the analysis two major categories of laminates are considered: 
(1) laminates with a constant nano tube volume fraction of 1% in which nano tube length 
is varied from 3 nm to 450 nm, and, (2) laminates with consisting of a constant a 
nanotube length of 400 nm in which nanotube volume fraction is varied from 1% to 20%. 
In the first category of laminates, shorter 3-nm nano tubes represent the current state-of- 
the-art, while longer 450-nm represent future technological goals. For the second 
category of laminates, volume fractions less than 5% represent the current state-of-the-art 



and volume fractions above 5%, while perhaps ideal, are in most cases too expensive or 
physically not practical. 

Hybrid Laminates 

Figure 3 shows a schematic of typical layups for uniform and hybrid angle-ply [+45/- 
45] s laminates. In the uniform laminate, the intrinsic molecular properties are constant 
throughout the laminate. By contrast, the hybrid laminate allows that the lamina 
properties may be alternated by varying the lamina molecular properties in a manner 
which adheres to laminate symmetry. Here, nanotube length, L, and volume fraction, v/, 
may vary about the midplane only if each layer located above the laminate midplane has 
an equivalent counterpart below the midplane. Table 3 shows a series of hybrid laminates 
of the [±45/0/90] ms family of laminates. 


Table 3. Q 

>uasi-isotropic hybrid [±45/0/901 ms laminate configurations. 

Laminate 

number 

Fixed lamina 
properties 

Laminate configuration 

1 

v/= 1% 

[(± 45 / 0 / 90)3 i= 1 0nm / (± 45 / 0 / 90 ) 3 L= % 

2 

v/= 1% 

[(± 45 / 0 / 90) 3 ' 400 nm /(± 45 / o / 90 ) 3 ' % 

3 

V/= 1% 

[(± 45 / 0 / 90 ) 3 L= x /(± 45 / 0 / 90) 3 i=400nm ] s 

4 

V/= 1% 

[(± 45 / 0 / 90 ) 3 L= x /(± 45 / 0 / 90) 3 ' 10nm ]s 

5 

L = 400 mn 

[(± 45 / 0 / 90)3 vf= 1 % /(± 45 / 0 / 90 ) 3 ,/= * ] s 

6 

L = 400 mn 

[(± 45 / 0 / 90)3 vf= 10 % /(± 45 / 0 / 90 ) 3 ,/= *] s 

7 

L = 400 mn 

[(± 45 / 0 / 90 ) 3 vf= * /(± 45 / 0 / 90) 3 vf= 10% ] s 

8 

L = 400 mn 

[(± 45 / 0 / 90 ) 3 vf= * / (± 45 / 0 / 90)3 vf= 1% ] s 


For each hybrid laminate configuration either nano tube length or volume fraction is 
fixed for the entire laminate. Also, the nanotube, length or volume fraction is fixed for 
groupings of lamina with a laminate denoted by the superscript within the laminate 
configuration. The nanotube length or volume fraction of the remaining groups of lamina 
may be assigned by varying the nanotube length and volume fraction parameters A and </>, 
respectively, thus allowing parametric studies. For example, the hybrid quasi-isotropic 
[ (±4 5 /0/9 0 ) 3 L 4<KI 11 m / (±45 /0/9 0 ) 3 L x ] s laminate, identified in Table 3 as laminate 2, has a 
constant volume fraction of 1% . The 3 outer (±45/0/90) groups on either side of the 
midplane are reinforced by 400-nm nanotubes, while the nanotube lengths of the inner 
(±45/0/90) groups may be assigned uniformly be varying the nanotube length parameter 
A. For the present study, A will be varied between 3 mn and 450 mn, and 4> will be varied 
between 1 % and 20 %. 



RESULTS AND DISCUSSION 


The buckling behavior of nanotube-reinforced polymer laminated composites will be 
assessed by presenting data showing the effects of stacking sequence and intrinsic 
molecular properties on the nondimensional parameters and buckling coefficients. 
Results are presented in this section for angle-ply [±0] ms ; quasi-isotropic [±45/0/90] ms ; 
axially-stiff [±45/02] ms ; ar| d transversely-stiff [±45/90 2 ] ms laminates. Also, to investigate 
the potential for optimizing the weight efficiency of such laminates, results depicting the 
effect of hybrid laminate configurations on the ratio of axial buckling load-to-nanotube 
weight percent are presented. 

Results showing the effects of varying nanotube length and volume fraction on the 
nondimensional bending orthotropy parameter f> for quasi-isotropic [±45/0/90] s laminates 
are presented in Figures 4(a) and 4(b). In Figure 4(a), nanotube length is varied uniformly 
for a constant 1% nanotube volume fraction. In Figure 4(b), nanotube volume fraction is 
varied uniformly where all plies are reinforced by 400-nm nanotubes. Both Figures 4(a) 
and 4(b) present results for quasi-isotropic laminates composed of aligned, randomly 
oriented, functionalized and non-functionalized nanotubes in various combinations. 

The results presented in Figure 4(a) indicate that for 1% volume fraction quasi- 
isotropic [±45/0/90] s laminates composed of functionalized and non-functionalized 
randomly-oriented nanotube-reinforced plies, the value of the nondimensional parameter 
j5 remains nearly constant over a range of nanotube lengths varying from 3 mn to 450 mn. 
In contrast, the data shows that for 1% nanotube volume fraction [±45/0/90] s laminates 
composed of aligned functionalized and non-functionalized nanotubes, the value of j5 
increases monotonically over the same range of nanotube lengths. The results for aligned 
non-functionalized CNT-reinforced laminates show a slightly larger value at the 
maximum nanotube length of 450 mn compared to the functionalized nanotube result. 

The results presented in Figure 4(b) indicate that for quasi-isotropic [±45/0/90] s 
laminates composed of 400-nm nanotubes, randomly-oriented functionalized and non- 
functionalized nanotube reinforcement yield a near constant response over a range of 
nanotube volume fractions increasing from 0% to 20%. By contrast, the results of Figure 
4(b) indicate that the [±45/0/90] s laminates composed of aligned functionalized and non- 
functionalized nanotubes show a nearly exponential increase in /3 as a function of volume 
fraction to a maximum at 5% nano tube volume fraction. As volume fraction increases 
beyond 5%, j3 remains nearly constant. 

The results of Figure 4, indicate the effect of intrinsic molecular properties on the 
orthotropy of nanotube-reinforced quasi-isotropic [±45/0/90] s laminate. In particular, 
results indicate in comparison to randomly-oriented nanotube-reinforced [±45/0/90] s 
laminates, aligned nanotube reinforcement yields larger values of the nondimensional 
bending orthotropy parameter j3. For both randomly-oriented and aligned nano tube 
reinforcement, functionalization has little effect on the degree of orthotropy. 

Figure 5 presents results that illustrate the effect on the nondimensional bending 
orthotropy parameter [:> due to varying nanotube length and volume fraction in axially- 
stiff [±45/02 ] s and [±45/02]6 S ; transversely-stiff [±45/902] s and [±45/90 2 ]6 S ; quasi- 
isotropic [±45/0/90] s and [±45/0/90] 6s ; and angle-ply [±45] s , [±45]6 S , [±30]6 S , and [±60]6 S 
laminates. The data shown in Figure 5(a) results from varying nanotube length for 
uniform laminates with 1% nanotube volume fraction. Figure 5(b) shows results of 



varying nanotube volume fraction for uniform laminates composed of 400-nm nanotubes. 
The results presented are for aligned, non-functionalized nanotubes. 

In Figure 5(a), the data shown indicates that for all laminates considered the 
nondimensional parameter /I increases monotonically as nanotube length in the uniform 
laminate increases from 3 mn to 450 mn. For the 3 mn nanotube length, the quasi- 
isotropic axially-stiff and transversely-stiff laminates share the maximum value of j> and 
the [±45] s and [±45]6 S angle-ply laminates share the minimum value. As the nanotube 
length increases uniformly to 100 mn, the values of the orthotropy parameter (5 for all of 
the laminates coalesce around unity. As the nanotube length approaches 450 mn, the 
maximum and minimum values of j> correspond to angle-ply [±45] family of laminates 
and quasi-isotropic [±45/0/90] 6s , respectively. 

The results shown in Figure 5(b) indicate a similar trend in the data for the laminate 
configurations, where a uniform nanotube length of 400 mn is maintained in each layer as 
nanotube volume fraction is varied from 0% to 20%. Here, the data coalesces around 
unity at about 1% volume fraction. At a nanotube volume fraction of about 5%, the 
maximum and minimum values of j3 correspond to angle-ply [±45] s and [±45]6 S laminates 
and [±45/0/90]6s, respectively. The values of j3 remains nearly constant over the range of 
volume fraction values from 5% to 20%. 

The results in Figure 5 indicate that for aligned, non-functionalized nano tube 
reinforcement the orthtropy for selected uniform laminates increases monotonically as 
nanotube length is increased for uniform 1% volume fraction laminates and as nanotube 
volume fraction is increased for uniform laminates composed of 400-nm nanotubes. The 
results show five distinct groupings of laminates that share orthotropy values: (1) [±45] s 
and [±45] 6s ; (2) [±30] 6s , and [±60] 6s ; (3) [±45/90 2 ] s , [±45/0 2 ] s , and [±45/0/90] s ; (4) 
[±45/0 2 ]6s and [±45/90 2 ]6 S ; and (5) [±45/0/90] 6s . These groupings of laminates are each 
represented by a line in the plots of Figures 5(a) and 5(b). The data in both Figures show 
that varying the intrinsic molecular properties significantly alters the degree of orthotropy 
among the unifonn laminates studied. Namely, trends in orthotropy associated with 
laminate configuration reverse themselves as nanotube length or volume fraction is 
increased. Thus, the maximum value of /! among the uniform laminates considered can 
be tailored by varying the intrinsic molecular properties of the nanotube used to reinforce 
the laminates. 

Nondimensional buckling coefficients for axial compression and shear loadings, K x 
and K s , respectively, for a quasi-isotropic [±45/0/90] 6s unifonn laminate are shown in 
Figure 6. Results are shown for both simply-supported and clamped boundary conditions. 
The nondimensional buckling coefficients are presented as functions of nanotube length 
for a constant 1% volume fraction for aligned, non-functionalized carbon nanotubes in 
Figure 6(a). While, Figure 6(b) shows nondimensional buckling results as a function of 
nanotube volume fraction for 400-nm aligned, non-functionalized CNT reinforcement. 

In Figure 6(a), the axial buckling coefficient and shear buckling coefficient, K x and 
K s , respectively, for the uniform quasi-isotropic [±45/0/90] 6s laminate both show a slight 
increase over the range of nanotube length from 3nm to 450 mn. The clamped and 
simply-supported results show the same trend. For both buckling coefficients K x and K s , 
the simply-supported results, represented by the dashed lines, have greater values than the 
clamped results. Values of the shear buckling coefficient K s are shown to exceed values 
of the axial buckling coefficient K x over the entire range of nanotube lengths. 



A similar trend is observed in the buckling results presented as a function of nanotube 
volume fraction shown in Figure 6(b). Here, the response of a uniform quasi-isotropic 
[±45/0/90]6 S laminate is considered for a 400-nm aligned, non-functionalized nanotube 
reinforcement by varying nanotube volume fraction from 0% to 20%. Both buckling 
coefficients K x and K s increase slightly over the range of nanotube volume fractions. The 
results for both the clamped and simply-supported boundary conditions show the same 
trends as functions of volume fraction. The simply-supported results, represented by the 
dashed lines, are greater than the buckling coefficient values for the clamped boundary 
conditions. Also, the value of the shear buckling coefficient K s is shown to exceed the 
values of the axial buckling coefficient K x for entire range of nanotube volume fractions. 

Buckling results presented in Figure 6 for the uniform quasi-isotropic [±45/0/90] 6s for 
aligned, non-functionalized CNT show that buckling response is a function of loading 
and boundary conditions. For a particular laminate, trends in the buckling results for axial 
compression and shear loadings are similar. Furthermore results for both simply- 
supported of clamped boundary conditions are similar. 

Buckling results for simply-supported uniform axially-stiff [±45/02] 6 S ; transversely- 
stiff [±45/90 2 ] 6s ; quasi-isotropic [±45/0/90] 6s ; and angle-ply [±45]6 S laminates reinforced 
with aligned, non-functionalized CNT are presented in Figures 7 and 8. Results for the 
axial buckling coefficient K x are presented in Figure 7 as a function of nano tube length 
for a constant 1% volume fraction. In Figure 8, the shear buckling coefficient results as a 
function of nanotube volume concentration are presented for a 400-nm CNT. 

Figure 7 shows that for 1% nanotube volume fraction, at a nanotube length of 3nm, 
the uniform angle-ply [±45]6 S laminate has the lowest axial buckling coefficient, while 
the uniform axially-stiff [±45/0 2 ]6 S , transversely-stiff [±45/90 2 ]6 S , and quasi-isotropic 
[±45/0/90]6 S laminates share the same value. As nanotube length increases the value of 
the buckling coefficients increase monotonically. At around 150 mn, the value of the 
buckling coefficient for the [±45]6 S laminate becomes largest. Finally, at 450 mn, the 
[±45]6 s laminate has the maximum value, the [±45/0 2 ]6 S and [±45/90 2 ]6 S share the same 
value, and the [±45/0/90] 6s has the lowest value. 

A similar trend for the shear buckling coefficient for 400-nm aligned, non- 
functionalized nanotube reinforcements over a range of nanotube volume fractions from 
0% to 20% is shown in Figure 8. At 0%, the [±45]6 S laminates have the shear buckling 
coefficient, while the uniform [±45/0 2 ]6 S , [±45/90 2 ]6 S , and [±45/0/90] 6s laminates share 
the same greater value. As nanotube volume increases the value of the shear buckling 
coefficients increase monotonically. At around 2% nanotube volume fraction, the value 
of the shear buckling coefficient for the [±45]6 S laminate becomes largest. Over the range 
of nanotube volume fractions from 5% to 20%, the [±45]6 S laminate has the maximum 
value, the [±45/0 2 ]6 S and [±45/90 2 ] 6s share the same value, and the [±45/0/90] 6s has the 
lowest value. 

The results of Figures 7 and 8 show that the intrinsic molecular properties may have a 
significant effect on the buckling response of certain nanotube-reinforced laminates 
exposed to axial compression and shear loadings. In particular, buckling coefficient 
values of the uniform angle-ply [±45]6 S laminate ranged from minimum to maximum in 
both axial compression and shear buckling response over a range of nanotube length and 
volume fraction values for aligned, non-functionalized carbon nanotube reinforcement. 
The increase in buckling load due to variation in nanotube proportion and size suggest the 



potential to tailor the buckling response of the laminate to various loading conditions 
using the molecular properties of the nanotubes. 

Figure 8 presents buckling results for the hybrid quasi-isotropic [±45/0/90] 6s family of 
CNT-reinforced laminates subjected to axial compressive loadings. These laminates are 
considered hybrid due to the fact that the molecular properties of the aligned, non- 
functionalized CNT reinforcement are varied through-the-thickness of the laminates in a 
manner similar to that shown in Figure 3. The laminates considered are listed in Table 3. 
Here, the buckling results are presented as a ratio of axial buckling load-to-nanotube 
weight percent. The data is normalized by the ratio of axial buckling load-to-nanotube 
weight percent ratio for a uniform quasi-isotropic [±45/0/90]6 S laminate. 

Figure 8(a) shows the ratio of buckling load-to-nanotube weight percent for hybrid 
laminates 1-4 in Table 3, where the nano tube length parameter A is varied from 3 mn to 
450 mn for a constant nanotube volume fraction of 1%. Here, the results are presented as 
a function of nanotube length. The data is nonnalized by a uniform quasi-isotropic 
[±45/0/90]6 S laminate with a constant 1% volume. As the nanotube length is increased, 
the hybridized layers in the laminate remain at a constant value denoted by the 
superscripts in the laminate configuration notation as described previously. The results in 
Figure 8(a) show that, in comparison to the uniform [±45/0/90] 6s laminate, hybrid 
laminates 2 and 3, where the hybrid layers are fixed at 400 mn, exhibit a higher buckling 
load-to-nanotube weight percent ratio value for a A = 10 mn. Laminates 1 and 4, where 
the hybrid layers are fixed at 10 mn, share buckling coefficient values with the uniform 
laminate at A = 10 mn. As the nanotube length in the variable layers increases, the 
normalized buckling load-to-nanotube weight percent ratio for hybrid laminates 2 and 3 
approaches unity, while the values for laminates 1 and 4 decrease to values less than 
unity. At the maximum nanotube length analyzed, 450 mn, laminate 1 has the the 
smallest buckling load-to-nanotube weight percent ratio. 

All laminates in Figure 8(a) have a constant nanotube volume fraction of 1%. 
Therefore, in laminates 2 and 3 the increase in axial buckling capacity relative to the 
volume of nanotubes achieved at A = 10 mn constitutes an enhancement in buckling 
performance without adding mass. Furthermore, utilizing hybrid layers away from the 
midplane appears to offer the greatest advantage. At A = 10 mn, laminate 2 shows the 
largest increase in the nonnalized axial buckling load-to-nanotube weight percent ratio. 

Figure 8(b) shows the normalized ratio of buckling load-to-nanotube weight percent 
for hybrid laminates 5-8 in Table 3, where the nanotube volume fraction parameter (p is 
varied from 1% and 20% for constant nanotube length of 400 mn. Here, the results are 
presented as a function of nanotube volume fraction. The data is normalized by a 
uniform [±45/0/90]6 S laminate composed of 400 mn nanotubes. The results for laminates 
5 and 8 coincide with unity at 1% volume fraction. The results for laminates 6 and 7 
coincide with unity at 10% volume fraction. This indicates that these laminates have the 
same volume fraction and buckling response as the uniform laminate at these points. 
However, for laminates 5 and 8, as the volume fraction parameter is increased, the 
normalized ratio of buckling load to weight percent nanotubes increases over the range of 
volume fraction values. At </> = 5%, laminate 5 has the maximum value at 1.7. Laminate 8 
reaches a value just below 1.7 at ^ = 5%. This slight difference is due to the fact that the 
fixed 1% volume fraction layers are placed at the midplane in laminate 8 and at the outer 
surface in laminate 5. More importantly the results shown reflects a 70% improvement in 



the buckling load-to-weight percent ratio due to the reduction in overall nanotube volume 
fraction in the laminate. For </> > 5% there is little change in the value of the buckling 
load-to-nanotube weight percent ratio. Also, for laminates 6 and 7, the load to weight 
percent ratio ranges from 0.2 at (p = 1% to 1.25 at </> = 20%. The maximum value for 
laminates 6 and 7 represents a 25% increase in the load to weight percent ratio. 

The results of Figure 8(b) demonstrate the fact that the axial buckling load-to-weight 
percent ratio can be improved by varying the nanotube volume fraction through-the- 
thickness of a laminate. Fixing the hybrid layers at 1% volume fraction yielded increases 
in the buckling load-to-nanotube weight percent ratio across a range of volume fractions 
from 1% to 20%. Fixing the volume fraction in hybrid layers at 10% yielded increases in 
the load to weight percent ratio across a range of volume fractions from 10% to 20%. 
Placing the hybrid layers away from the midplane offered a slight improvement in the 
case of the 1% volume fraction hybrid layers, but not so with the 10% hybrid layers. The 
hybrid laminates with fixed 1% nanotube volume fraction layers offered the maximum 
buckling load-to-nanotube weight percent improvement of 70% compared to the uniform 
laminate. The results suggest that the weight efficiency of the nanotube-reinforced 
laminates can be improved by optimizing the nanotube volume fraction with respect to 
laminate configuration. Furthermore, the optimal buckling load-to-nanotube weight 
percent ratio can be achieved using low nanotube volume fractions. 


SUMMARY AND CONCLUSIONS 

An analysis to predict the effect of intrinsic molecular properties on the buckling 
response of polyethylene composites reinforced by CNT subjected to axial compression 
and in-plane shear loadings has been conducted. The results indicate the potential for 
tailoring the macromechanical response of the CNT-reinforced laminates built by 
layering of plies or lamina of polyethylene matrix material reinforced by CNT suspended 
in the matrix. The buckling response of layered polyethylene CNT-reinforced laminates 
was demonstrated by developing a model to predict the critical buckling response of the 
symmetric laminates based upon the intrinsic molecular properties of the nanotube 
length, volume fraction, orientation, and functionalization. 

Results indicate in comparison to uniform randomly-oriented CNT-reinforced quasi- 
isotropic [±45/0/90] s laminates, aligned nanotube reinforcement yields larger values of 
the nondimensional bending ortho tropy parameter j3. For both randomly-oriented and 
aligned nanotube reinforcement functionalization has little effect on the degree of 
orthotropy. For aligned, non-functionalized nanotube reinforcement the orthtropy for 
selected uniform laminates increases monotonically as nanotube length is increased for 
1% volume fraction laminates and as nanotube volume fraction is increased for laminates 
composed of 400-nm nanotubes. The results show distinct groupings of laminates that 
share orthotropy values. Varying the intrinsic molecular properties significantly alters the 
degree of orthotropy among the uniform laminates studied. Namely, trends in orthotropy 
associated with laminate configuration reverse themselves as nanotube length or volume 
fraction is increased. Increasing volume fraction of shorter nanotubes or increasing 
nanotube length for low volume fraction composites gives rise to trends in laminate 
bending orthotropy similar to those encountered in fiber-reinforced polymer composites. 



Buckling results show that buckling response is a function of loading and boundary 
conditions. For a particular laminate, trends in the buckling results for axial compression 
and shear loadings are similar. Furthennore results for both simply-supported of clamped 
boundary conditions are similar. Results show that the intrinsic molecular properties may 
have a significant effect on the buckling response of certain CNT-reinforced laminates. In 
particular, among a set of selected laminates configurations, buckling coefficient values 
of the unifonn angle-ply [±45]6 S laminate ranged from minimum to maximum in both 
axial compression and shear buckling response over a range of nanotube length and 
volume fraction values for aligned, non-functionalized CNT reinforcement. The increase 
in buckling load due to variation in nanotube proportion and size suggest the potential to 
tailor the buckling response of the laminate using the molecular properties of the 
nanotubes. 

Finally, analysis of the hybrid laminate configurations showed increases in axial 
buckling load may be achieved without adding mass by varying the nanotube length 
through- thickness. Futhermore, utilizing hybrid layers away from the midplane appears to 
offer the greatest advantage. For A= 10 mn, [(±45/0/90) 3 £ “ 400nm /(±45/0/90)3 £-x ] s shows a 
15% increase in the normalized buckling load-to-nanotube weight percent ratio compared 
to a unifonn [±45/0/90] 6s laminate. By contrast, [(±45/0/90)/ _x /(±45/0/90)3 / ' 400nni ] s 
shows only a 7% increase. Results indicate that while longer nanotubes may be optimal 
for mechanical properties in a generic sense, a mix of long and short nanotubes may be 
optimal for enhancing buckling perfonnance. Other results indicate that the axial 
buckling load-to-nanotube weight percent ratio can be improved by hybridizing the 
nanotube volume fraction of the laminate. The hybrid [(±45/0/90) 3 1 ^ 1 /o /(±45/0/90) 3 % 

laminate offered the maximum buckling to weight percent improvement of 70% at <j> = 
5% compared to the uniform quasi-isotropic laminate indicating that the volume fraction, 
and therefore cost, maybe optimized with respect to laminate configuration to enhance 
buckling perfonnance at nano tube volume fractions less than 5%. 
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FIGURES 



Note: For k = 1,7V, Zk-i - = constant 

Figure 1. Geometry of //-layered laminate. 



Figure 2. Geometry of infinitely long symmetric laminate subjected to axial 
compression and shear loadings. 
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Note: Vf, nanotube volume fraction; 0, nanotube orientation angle; and L, nanotube length 
Figure 3. Schematic of typical [±45] s uniform and hybrid symmetric laminates. 
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(a) 1% nanotube volume fraction 



(b) 400-nm nanotube length 


Figure 4. Effect of intrinsic molecular properties on non-dimensional bending ortho tropy 
parameter, (3, for [±45/0/90] s laminates. 






(a) 1% nanotube volume fraction 



Figure 5. Effect of aligned, non-functionalized nano tube length and volume fraction on the 
non-dimensional bending orthotropy parameter, (3, for selected laminates. 
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(b) 400-nm nanotube length 

Figure 6. Effect of aligned, non-functionalized nano tube length and volume fraction on the 
non-dimensional buckling coefficients, K x and K s , for [±45/0/90] 6s laminates subjected to 
axial compression and shear loadings. 






(a) 1% nanotube volume fraction 

Figure 7. Effect of aligned, non-functionalized nano tube length on the non-dimensional axial 
buckling coefficient K x for selected laminates subjected to axial compression loadings with 
simply-supported boundary conditions. 


. _L _j„ 



volume fraction, % 

(b) 400-nm nanotube length 

Figure 8. Effect of aligned, non-functionalized nano tube volume fraction on the non- 
dimensional shear buckling coefficient K s for selected laminates subjected to shear loadings 
with simply-supported boundary conditions. 







• Laminate 1: [(±45/0/90) 3 L=10nm /(±45/0/90) 3 L=A ]s 

♦ Laminate 2: [(±45/0/90) 3 L=4 °° nm /(±45/0/90) 3 L=A ] s 
a Laminate 3: [(±45/0/90) 3 L=A /(±45/0/90) 3 L=400nm ] s 
■ Laminate 4: [(±45/0/90) 3 L=A /(±45/0/90) 3 L=10nm ] s 
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Figure 9. Effect of hybrid laminate construction on the ratio of axial buckling load-to- 
nanotube weight percent for aligned, non-functionalized quasi-isotropic laminates subjected 
to axial compression loadings with simply-supported boundary conditions. 



